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Neutron spin-flip scattering observations of a resonance in the superconducting state is often taken 
as evidence of an unconventional superconducting state in which the gap changes sign A(k + Q) = 
— A(fe) for momentum transfers Q which play an important role in the pairing. Recently questions 
regarding this identification for the Fe-pnictide superconductors have been raised and it has been 
suggested that A(k + Q) = A(k). Here we propose that inelastic neutron or x-ray scattering 
measurements of the spectral weight of a phonon of momentum Q can distinguish between these 
two pairing scenarios. 



Determining the superconducting gap structure in the 
Fe-based superconductors has proven difficult due to 
the complicated multishcctcd, multiorbital nature of the 
Fermi surface, as well as problems with fabricating clean 
films and interfaces. At present the most popular can- 
didate for pairing in most of these materials is the 
s~! state which changes sign between hole and elec- 
tron pockets and arises naturally from a spin-fluctuation 
interaction[l-3]. However, there are few phase-sensitive 
tests of the relative sign change of Ak between the dif- 
ferent Fermi surfaces. Josephson tunneling experiments 
[4, 5], Fourier transform STM in a magnetic field [6], and 
the observation of a resonance mode in inelastic neutron 
scattering experiments[7-10] have been offered as limited 
evidence for such a sign change. 

Recently, an electron-phonon interaction enhanced by 
orbital fluctuations has been proposed as the mechanism 
which is responsible for pairing in the Fe-based supercon- 
ductors [11-13]. This interaction leads to a so-called s ++ 
gap which is positive on both the hole and electron Fermi 
surfaces. While in both cases, the gap can be anisotropic 
and accidental nodes can even occur, the relative sign 
of the gap on regions of the hole and electron Fermi sur- 
faces which have the same dominant orbital weight will be 
negative if the pairing is driven by spin-fluctuations and 
positive if driven by phonon-orbital-fluctuations. Thus it 
is important to determine the relative sign between A(fc) 
and A(k + Q) on the hole and electron Fermi surfaces, 
respectively, for momentum transfers Q which connect 
such regions. 

One test of this is the occurrence of a spin resonance 
peak in the neutron spin-flip scattering in the supercon- 
ducting state [14]. The BCS coherence factor which en- 
ters the spin-flip scattering is 

If A(fc + Q)A(fc) \ 

2 ^ E(k + Q)E(k) J' 1 ' 



Here E(k) = A 2 (fc) is the quasi-particle energy 

and near threshold, where goes to zero, this coher- 
ence factor goes to 1 if A(fc + Q) and A(fc) have op- 
posite signs and zero if they have the same sign. A 
resonance peak is expected in the former case [15, 16]. 
However, arguments [17] have been given that for the 
case of BaFei.85Coo.i5As2 the experimentally observed 
peak [7-10] may be consistent with an s ++ gap, if the 
scattering rate collapses sufficiently rapidly in the super- 
conducting state. Thus one would like to find additional 
experimental probes which provide a way of determining 
the relative sign of the gap between the hole and electron 
pockets. 

As is known, there exist anomalies in the phonon spec- 
tral weight of conventional superconductors. Various au- 
thors [18-20] have shown that if the phonon frequency £l q 
is greater than twice the gap A, but not too far above, 
significant spectral weight is transferred below 2 A into a 
resonance. In addition, the phonon peak shifts to higher 
frequencies and broadens. These effects have been con- 
vincingly observed in neutron scattering experiments on 
the borocarbide superconductor YNi 2 B 2 C [21]. This be- 
havior reflects the fact that in the superconducting state, 
the coherence factor that enters the coupling of a phonon 
to electron polarization fluctuations is similar to Eq. (1) 
except that it has a plus sign. In this case the coher- 
ence factor will go to 1 at threshold if A(fc + Q) and 
A(fc) have the same sign. This will lead to an anoma- 
lous resonance and hardening of the phonon peak for an 
s ++ gap. Here we point out that for an s" 1 gap with 
A(fc) = — A(k + Q), the anomalous resonance is absent 
and the phonon peak shifts down (softens) slightly. Thus 
neutron and x-ray scattering [22] measurements of the 
phonon spectral weight at specific Q values can provide 
information on the relative sign of the gap [23] . 

To explore the use of this effect to determine the rel- 
ative gap signs in the Fe-based materials, we will use 
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the same 5-orbital tight binding fit of the LDA band- 
structure that we previously used to discuss the possible 
occurrence of a neutron scattering spin resonance [24]. 
In this case, with an electron doping x = 0.125 one finds 
the Fermi surfaces and orbital occupations illustrated in 
Fig. 1. Here there are two hole Fermi surfaces a\ and 
oli around the T point and two electron Fermi surfaces 
Pi and /?2 around the X(ir,0) and Y(0,tt) points of the 
unfolded (lFe/cell) Brillouin zone. We will use an or- 
bital notation £ = 1, 2, 3, 4, 5 to denote the d xz , d yz , d xy , 
d x 2_ y 2 and d 3z 2_ r 2 orbits, respectively, where the axes 
are along the Fe-Fe bonds. The dominant orbital weight 
of the Bloch states on the Fermi surfaces is indicated in 
Fig. 1. 




FIG. 1. (Color online) Fermi surfaces for a five-orbital model 
of the Fe-pnictide superconductors. The colors indicate the 
dominant orbital character (red=d xz , green=d yz , b\ue=d xy ). 
Here Qi connects d xz parts of the hole 012 and d yz parts of 
the electron /3i Fermi surfaces and Q2 connects d xz parts of 
ct\ to the d xy parts of j3\. 



As Kontani and Onari [11] have noted, local Fe-ion 
displacements couple to orbital fluctuations. The polar- 
ization susceptibility associated with these orbital fluc- 
tuations enters the phonon self-energy and is reflected in 
the phonon spectral weight. Here we examine the change 
in the phonon spectral weight — ^Im D(Q,oj) associated 
with the Qi = (7r,0.097r) and Q2 — (0.86-7T, 0) momenta 
shown in Fig. 1 when the system goes from a normal 
state to a superconducting state with a nominal s ++ gap 
or an s" 1 gap. We begin with Qi & n d will assume that 
there is an electron-phonon coupling 1712 for scattering 
an electron between the dominantly d xz -\ike states on 
the «2 hole Fermi surface to the d yz -like states on the j3% 
electron Fermi surfaces. In the following, the real part 
of the normal state phonon self-energy will be absorbed 
to give the dressed normal state phonon frequency CIq 1 . 
Then at low energies, the imaginary part of the phonon 



self-energy for momentum Qi dominantly comes from 
scattering between the nearly nested d xz -lihe states on 
«2 and the d yz -like states on j3\. Similarly, the change 
between the superconducting and normal state phonon 
self-energy is also determined by the d xz -d yz orbital fluc- 
tuations. In this case the important orbital contribution 
to the phonon self-energy is 

n(g )W ) = -|. 9l2 | 2 (xSkfow) +x?Mq,")) (2) 

with 

X1212M = fx°(9,w) [1 + U c x°(q^)V 1 } (3) 
I J 1212 

and similarly for Xam- Here the Coulomb interaction 
matrix is 



U l x = 4 = h = h 

-U' + 2 J tx = 4 ± I2 = ii 

j ex = £ 2 ^i 3 = h 

J' ti = U±i 2 = h 



(4) 



and 



x {G^fc + q)G v (k) - F»(-k - q)F»{k)}. 



(5) 



The notation follows Ref. [25] with U and U' the intra- 
and inter-orbital Coulomb interactions, J the exchange 
coupling and J' the pair hopping interactions. We will 
consider a typically rotationally invariant set of inter- 
action parameters which are similar to those previously 
used to study the occurrence of a spin resonance in the 
s + ~ superconducting state [24]. These are such that 
U = 4 in units of the hopping tj 1 , J = £7/16, U' = U-2 J, 
and J' — J [26] . In the superconducting state the normal 
and anomalous Gor'kov Green's functions are given as 



G"(jfc) 



F"(k) = 



A M (fc) 



(6) 



with E n 



For the normal state calcu- 



lations, the gap is set to zero in Eq. (5). 

Here, because we are calculating the orbital- 
fluctuations, the sign between the normal and anoma- 
lous Green's function products in Eq. (5) is negative. It 
would be positive if we were studying the spin suscepti- 
bility and in addition the orbital interaction matrix U c 
would be changed to — II s . In Eq. (5) we use a nota- 
tion in which k = (k, w n ) and q = (q, (u m — > u) + id)). 
The calculations were carried out for a low temperature 
T = 0.05 in units of t\ x . 

Figure 2 shows the imaginary part of X^212(Qi Cl ') versus 
u with q — Qi for the normal state (solid), the supercon- 
ducting states with an s ++ (A Ql = A„ 2 = = A i a 2 = 
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FIG. 2. (Color online) The imaginary part of Xi2i2(<3i 5 w) 
versus to for the normal state (red, solid), the s ++ supercon- 
ducting state with Ao(fc + Qi) = Ao(fc) = Ao (blue, dashed) 
and the s H state with Ao(fc + Qi) = — Ao(fc) = Ao (green, 
dot-dashed) . 



A ) gap (dashed) and an s H (A Ql = A Q2 = — A^ = 
— Ap 2 = A ) gap (dash-dot) [27]. The imaginary part 
of the orbital response is similar to what was found for 
the spin susceptibility with the s H and s ++ gaps inter- 
changed. Thus, if the phonon-orbital-fluctuations pro- 
vide the dominant pairing interaction in the Fe-based 
superconductors then the gap will have s ++ character 
and one will find an anomalous resonance in the orbital 
susceptibility in the superconducting state rather than in 
the spin susceptibility. 

Although it would be interesting to have a direct scat- 
tering probe of the orbital susceptibility spectral weight 
shown in Fig. 2, one can still access its structure by mea- 
suring the phonon spectral weight. Using the orbital sus- 
ceptibility, we can evaluate the phonon self-energies that 
enter the normal and superconducting phonon progaga- 
tors [28], 



and 



Ds 1 (q,uj) = 



2il t 



2n c 



- Rc5~tt{Q 



^-iImniv(Qi,w) 



(7) 



l, U) - iImIls(Qi,uj). 

(8) 

Here <5IT = II5 — ITat and lis and II n are the supercon- 
ducting and normal phonon self-energies obtained from 
Eq. (2), with the gap set to zero in Eq. (5) for the 
normal case and set to A(fc) = A(fc + Q) — An and 



A(fc) = -A(k + Q) = A for the s 
ducting states, respectively. 



++ 



and s + supercon- 



o 

a 

'3 

a 
a 



(a) 



normal 



1 - s T 
- s + 



a 



a 
a 



(b) 



Re8n s .*/£2 Ql 
Re5n s .-/£2 Ql 




O)/A 



FIG. 3. (Color online) (a) The imaginary part of the phonon 
damping n(Qi,a;) versus cj for the normal state (red, solid), 
the s ++ (blue, dashed) and s + ~ (green, dot-dashed) super- 
conducting gaps, (b) The real part 8U(Qi,uu) (see text) for 
the s ++ (red, solid) and s" 1 (blue, dashed) superconducting 
gaps. 

The phonon damping — lmIl(Q, u))/Qq 1 versus w/An 
is plotted in Fig. 3a. Here we have taken — 2.5A 
and set 2\g 12 \ 2 N(0)/£lQ 1 = 0.2. As CIq 1 approaches 2A , 
the effects we are discussing are enhanced and diminish 
when £Iq becomes large compared with 2Arj. As seen 
in Fig. 3a, for an s ++ gap there is an enhanced phonon 
damping above 2Arj. This arises from the increase in the 
density of states in this region due to the opening of the 
superconducting gap. This enhancement also depends 
upon the coherence factor and as seen for the case, 
if the coherence factor vanishes, the onset of the damping 
rises slowly as uj exceeds 2An and is smaller than in the 
normal state. 

In Fig. 3b, we have plotted the difference between 
the real part of the phonon self-energy in the supercon- 
ducting and the normal states asociated with the xz- 
yz scattering processes. The s ++ case shows the famil- 
iar hardening of the phonon that occurs when £Iq 1 ex- 
ceeds 2Arj. That is, in the s ++ superconducting state, 
Re(5n(Qi, uS) is positive for to > 2Aq. As we will see when 
we look at the phonon spectral weight, the negative dip 
for u> < 2 An leads to what Allen et al.[20] have called a 
"vibrational/superelectronic" resonance. The structure 
of Re<5n(Qi, cj) is quite difference in the state. Here 
for cj > 2A , SReU(Qi,uj) is negative and this leads to a 
softening of the ^Iq 1 phonon. In addition, there will not 
be a resonance feature below 2Aq in the phonon spectral 
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weight. 

Using the results for the self-energy shown in 
Fig. 3, we have plotted the phonon spectral weight 
— ^ i ImD((5i, uj) versus u>/A in Fig. 4a. The solid 
line shows the result for the normal state while the long 
dashed and the dash-dot curves are for the superconduct- 
ing s ++ and s H states, respectively. Here one sees that 
for an s ++ gap there is the expected resonance just below 
2Ao and the phonon peak has shifted up in frequency and 
broadened. For the s H gap, the anomalous resonance 
below 2Ao is absent and the phonon peak has slightly 
softened as well as narrowed. Similar results for the Q2 
phonon spectral weight are shown in Fig. 4b. 
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— ^-lmD(Q, uj) versus lo/Ao for wavevectors (a) Qi and (b) 
Q2 for the normal state (red, solid), the s ++ (blue, dashed) 
and s" 1 (green, dot-dashed) superconducting gaps. Note that 
we have added a finite damping in the calculation of the or- 
bital susceptibility x° Tb (<lyUj). 



In RPA calculations [3, 25, 29] there were parameter 
ranges in which the -like gap was found to have nodes 
on the (3 Fermi surfaces. In this case, the nominally s H 
gap will have the same sign on the parts of the Fermi 
surfaces connected by Q2- However, they would differ for 
Q\. Thus, it will be important to measure the phonon 
spectral weight for various near nesting Q values. How- 
ever, one will need phonon modes that couple to the rel- 
evant orbital scattering and in addition have frequencies 
Qq of order the sum of the magnitudes of the super- 
conducting gaps |A(fc + Q)\ + |A(fc)| associated with the 
nearby nested regions. In principle, with sufficient q and 
uj resolution and sensitivity, neutron and x-ray scattering 
offer the possibility of mapping out the relative structure 



of the gap between various parts of the Fermi surfaces in 
Fe-based as well as other unconventional superconduc- 
tors. 
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